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Diamagnetic Suppression of Component Magnetic 
Reconnection at the Magnetopause 

M. Swisdak, 1 B. N. Rogers, 2 J. F. Drake, 1 and M. A. Shay 1 

Abstract. We present particle-in-cell simulations of collisionless magnetic reconnection 
in a system (like the magnetopause) with a large density asymmetry across the current 
layer. In the presence of an ambient component of the magnetic field perpendicular to 
the reconnection plane the gradient creates a diamagnetic drift that advects the X-line 
with the electron diamagnetic velocity. When the relative drift between the ions and elec- 
trons is of the order the Alfven speed the large scale outflows from the X-line necessary 
for fast reconnection cannot develop and the reconnection is suppressed. We discuss how 
these effects vary with both the plasma (3 and the shear angle of the reconnecting field 
and discuss observational evidence for diamagnetic stabilization at the magnetopause. 



1. Introduction 

Magnetic reconnection, the breaking and reforming of 
magnetic field lines with a concomitant transfer of energy 
from the field to the surrounding plasma, is thought to drive 
such phenomena as solar flares, magnetospheric substorms, 
and tokamak sawtooth crashes. Despite the suspected link- 
age, the resistive magnetohydrodynamic (MHD) model of 
reconnection due to Sweet [1958] and Parker [1957] (the fa- 
vored theoretical model for many years), predicts reconnec- 
tion rates in each of these systems that are typically orders 
of magnitude too low to be physically relevant. The dis- 
crepancy arises because the Sweet-Parker reconnection rate 
varies with the Spitzer resistivity which in turn depends on 
particle interactions that are rare in these collisionless plas- 
mas. 

Recent work suggests that terms important at small 
length scales but usually ordered away in resistive MHD, 
notably the Hall term in Ohm's law, lead to reconnection 
rates that are consistent with the observations [Shay et al., 
1999]. With the addition of this new physics the speed of 
the electron flows away from the reconnection site (the X- 
line) is no longer bounded by the ion Alfven speed, instead 
scaling inversely with the width of the current layer formed 
by the reconnecting magnetic fields. As a result, the out- 
ward electron flux remains large even as the width, con- 
trolled by non-ideal effects, becomes small. Downstream 
from the reconnection site the super- Alfvenic electrons slow 
to rejoin the ions and expand outwards in a wide layer remi- 
niscent of the Petschek [1964] model. This behavior has been 
confirmed by a variety of numerical simulations: two-fluid, 
hybrid (fluid electrons and particle ions), and full particle. 
In the GEM (Geospace Environment Modeling) reconnec- 
tion challenge each of several independent codes modeling 
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an identical (two-dimensional) system found features simi- 
lar to those described above [Birn et al, 2001 and references 
therein] . Laboratory experiments in the relevant regimes are 
challenging, but recent results [Brown, 1999; Ji et al., 1999] 
support some aspects of this picture. 

Although the GEM challenge established a possible mech- 
anism for fast collisionless reconnection, it did so for a rela- 
tively simple geometry. Experimental evidence suggests that 
in nature reconnection is likely to be more complex. For 
instance, at the Earth's dayside magnetopause the magne- 
tosphere, a region of low density but strong magnetic field, 
abuts the magnetosheath and its high density but weaker 
magnetic field. Still, early single-point spacecraft measure- 
ments of the magnetopause [Paschmann et al, 1979; Son- 
nerup et al., 1981] detected field signatures and particle dis- 
tribution functions suggesting the occurrence of reconnec- 
tion. A more recent study using data from multiple space- 
craft [Phan et al., 2000] found, at least in one case, the 
bidirectional jets that are a reconnection hallmark. Simi- 
lar cross-field pressure gradients occur in laboratory fusion 
experiments where their associated diamagnetic drifts are 
thought to control the onset of reconnection during the saw- 
tooth crash [Levinton et al., 1994]. 

The first theoretical explorations of asymmetric recon- 
nection [Levy et al., 1964; Petschek and Thome, 1967] were 
based on incompressible MHD and predict an outflow re- 
gion comprising a combination of slow MHD shocks and slow 
and intermediate MHD waves. Two-dimensional computer 
simulations using compressible MHD [Hoshino and Nishida, 
1983; Scholer, 1989], hybrid [Lin and Xie, 1997; Omidi et 
al, 1998; Krauss-Varban et al., 1999; Nakamura and Sc- 
holer, 2000] and particle [Okuda, 1993] codes have addressed 
these predictions. Common features include the preferential 
growth of the magnetic island (O-line) towards the magne- 
tosheath and a pronounced density drop on the magneto- 
spheric side of the boundary layer. Also, consistent with 
the picture of Levy et al. [1964], the current layer abutting 
the magnetosheath is generally stronger than that bordering 
the magnetosphere and takes the form of a series of discon- 
tinuities. There is disagreement, however, concerning the 
precise form of these discontinuities, possibly because the 
nonlocality associated with the motion of collisionless par- 
ticles parallel to the shock front makes the two-dimensional 
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(2-D) simulations more complicated than the 1-D models. 
Simulations have also been performed with an initial out-of- 
plane component of the magnetic field (nonzero B y in GSM 
coordinates) [Hoshino and Nishida, 1983; Krauss-Varban et 
al., 1999; Nakamura and Scholer, 2000]. This guide field 
tends to slow, but not completely suppress, reconnection 
and alter the structure of the shock transitions. 

In our investigations we begin with a simple magne- 
topause equilibrium model with an ambient pressure dif- 
ference and perform two-dimensional fully kinetic numerical 
simulations. We reproduce the qualitative features seen in 
previous work but also find that diamagnetic drifts produced 
by the pressure gradient advect the X-line. When the dia- 
magnetic velocity, v* = — (c/qnB 2 )Vp X B (c is the speed 
of light; q, n, and p are the species charge, density, and pres- 
sure; B is the magnetic field), is comparable to the Alfven 
speed, reconnection is completely suppressed. These dia- 
magnetic effects have not been previously reported in the 
context of magnetopause reconnection. They are ordered 
out of the MHD model [Hoshino and Nishida, 1983] and 
were not seen in the hybrid simulations, perhaps because the 
spatially localized resistivity used to trigger fast reconnec- 
tion effectively prohibits X-line advection. In spite of this, 
diamagnetic effects can be significant at the magnetopause. 
Taking n ~ 10 cm" 3 , B ~ 10~ 4 G, k B T ~ 100 eV, and 
gradients of order an inverse ion inertial length, one finds 
v*/va ~ 1. In the fusion context, previous analytic work 
has suggested that diamagnetic drifts can stabilize the tear- 
ing mode and lower the rate of reconnection [Biskamp, 1981; 
Rogers and Zakharov, 1995]. But since this work considered 
a reduced MHD model (equivalent to the limit of a large 
guide field) its applicability to magnetospheric applications 
has not been established. 

Unlike the configuration in the magnetotail, the magnetic 
field on opposite sides of the magnetopause is usually not 
equal in magnitude and anti-parallel. Reconnection in such 
a system does not occur at a unique spatial location, even 
in a model with a 1-D equilibrium, since components of the 
magnetic field reverse direction at any location across the 
current layer. Because of the possibility that reconnection 
can occur at multiple spatial locations, it has been suggested 
that the magnetopause magnetic field is stochastic [Galeev 
et al, 1986; Lee et at, 1993]. The exploration of reconnec- 
tion in a 2-D model must then be carried out with care since 
stability at a particular surface does not necessarily imply 
stability at all surfaces. In exploring the impact of diamag- 
netic drifts on reconnection we therefore must address both 
stability at a single plane and explore which plane is ex- 
pected to dominate the dynamics of a full 3-D system. We 
suggest on the basis of analytic arguments and simulations 
that the strongest reconnection occurs at the surface where 
the reversed field components have equal magnitudes. 

In section 2 of the paper we present our computational 
scheme and initial conditions. Section 3 is an overview of 
reconnection at the magnetopause, both without and with 
a diamagnetic drift. Section 4 discusses how varying the 
strength of the out-of-plane field changes the diamagnetic 
effects while section 5 addresses the question of determin- 
ing the dominant reconnection plane. Finally, in section 6 
we summarize our results and discuss their implications for 
understanding magnetic reconnection at the magnetopause. 

2. Computational Methods 

Although diamagnetic drifts are present in two-fluid and 
hybrid models, a complete calculation in the collisionless 
limit requires a careful, and assumption-filled, treatment of 



the full pressure tensor. To avoid this difficulty, we sacrifice 
the benefits of a fluid simulation for a more computationally 
demanding, but mathematically straightforward, full parti- 
cle description. 

2.1. The Code 

The simulations are done with p3d, a massively parallel 
kinetic code that can evolve up to ~ 10 9 particles on current 
Cray T3E and IBM SP systems [Zeiler et al, 2002]. The 
Lorentz equation of motion for each particle is evolved by a 
Boris algorithm (E accelerates for half a timestep, followed 
by a rotation of v by B, and then the other half-step acceler- 
ation by E) . The electromagnetic fields are advanced in time 
with an explicit trapezoidal-leapfrog method using second- 
order spatial derivatives. Poisson's equation constrains the 
system; if V • E 7^ 4-7rp, a multigrid algorithm recursively 
corrects the electric field. Although the code permits other 
choices, we work with fully periodic boundary conditions. 

The code is written in normalized units: lengths to the ion 
inertial length c/io pi — d», times to the inverse ion cyclotron 
frequency fij/, velocities to the Alfven speed va, masses to 
the ion mass rrii, and temperatures to m,iv\. Unless other- 
wise noted, the asymptotic reconnecting field is assumed to 
have a magnitude of 1 on both sides of the magnetopause 
(the reason for this choice will be made clear in section 5) 
and the asymptotic density in the magnetosheath is taken 
to be 1. 

To conserve computational resources, yet assure a suffi- 
cient separation of spatial and temporal scales, we take the 
electron mass to be 0.005 = 1/200 and the speed of light to 
be 20. Our simulations are performed in a box of 1024 x 1024 
gridpoints that is 25.6 on a side, corresponding to a grid scale 
of 0.025 and so ?»3 gridpoints per electron inertial length. 
A typical timestep is 1.5 x 10~ 3 . 

2.2. Initial Conditions 

Reconnection occurs in the x — y plane in our coordi- 
nate system. For reference, this is equivalent to the z — it- 
plane in GSM coordinates. The initial equilibrium is con- 
structed from two (modified) Harris current sheets centered 
at Ly/A and 3L y /4, where L y is the box size in the y direc- 
tion; to allow periodicity the current is parallel to z in one 
sheet and anti-parallel in the other. The current sheets pro- 
duce a magnetic field B x = tanh[(j/ — L y /4)/wo] — tanh[(jy — 
3L H /4)/wo] — 1, where wo = 0.5. The density profile is sim- 
ilar, n = no(tanh[(y — L y /4)/wo] — tanh[(y — 3L y /4)/wo]) + 
n m in, where, unless otherwise specified, no = 0.45 and 
n m in =0.1. This profile gives an asymptotic density of n — 1 
in the magnetosheath and 0.1 in the magnetosphere. The 
initial ion and electron temperatures are T% = 2 and T e = 1, 
implying thermal speeds of » 1.4 and 14, respectively. The 
guide (out-of-plane) field B z is assigned a specific asymp- 
totic value B z o on the magnetosheath side and is calculated 
elsewhere by assuming pressure balance. Note that at the 
reversal surface, where B x — 0, pressure balance implies 
that d y n(Ti + T e ) oc B z d y B z , and hence a density gradient 
must be accompanied by a gradient in B z . At t = a 5% 
perturbation in the magnetic flux function places the X- and 
O-line in each sheet. 

Unlike a conventional Harris sheet, this initial state is not 
a kinetic (Vlasov) equilibrium. Although it has been shown 
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[Quest and Coroniti, 1981] that the presence of a guide field 
reduces the growth rate of the tearing mode instability in the 
linear regime, our initial perturbation is large enough that 
the system begins to reconnect nonlinearly before any other 
perturbations become important. We must also emphasize 
that our system characterizes the magnetopause only in a 
rough approximation. In particular, magnetospheric plasma 
is typically not isothermal, but instead has temperature gra- 
dients comparable in magnitude, but anti-parallel, to those 
of the plasma density. If this effect is included the over- 
all pressure gradient, as well as the diamagnetic effects we 
discuss here, will be smaller. 

3. Simulation Overview 

In the simplest description of fast collisionless reconnec- 
tion three length scales are relevant: the system size and 
the ion and electron inertial lengths. (Depending on the 
strength of the initial out-of-plane field, the ion Larmor ra- 
dius can also enter as a parameter [Kleva and Drake, 1995].) 
The dynamic couplings between these scales, which must ex- 
ist for fast reconnection to occur, are discussed in detail by 
Biskamp et al. [1997] and Shay et al. [1998]. At large scales 
both species are frozen into the magnetic field and flow to- 
wards the X-line with a velocity Vi„ that is some fraction of 
the upstream Alfven speed. Roughly c/uu P i away from the 
current sheet the inertia term in the ion equation of mo- 
tion becomes important and the ions, but not the electrons, 
decouple from the field. Deflected outwards, the ions flow 
Alfvenically away from the X-line while the electrons travel 
further inward until they too decouple and travel outward 
at velocities near the electron Alfven speed. Downstream, 
the electrons slow to join the ionic outflow and reform an 




MHD fluid. In what follows we show that diamagnetic drifts 
can disrupt this structure by strongly altering the MHD out- 
flows, thereby inhibiting the coupling of the ions to the X- 
line and halting the formation of a large-scale reconnection 
geometry. 

3.1. Anti-parallel Reconnection 

To document certain features and establish contact with 
previous work, Figure 1 shows the X-line structure from 
a magnetopause simulation similar to those described by 
Krauss-Varban et al. [1999] and Nakamura and Scholer 
[2000]. The upper section of each panel {y/di > 3.2) con- 
tains low density plasma and a large magnetic field like the 
magnetosphere, while the lower section has a high density 
plasma and a weak field like the magnetosheath. There are 
a few differences between the initial conditions for this sim- 
ulation and the case with a guide field discussed in section 
2.2 (and used in the rest of this paper). The reversed field 
is given by a shifted hyperbolic tangent such that the mag- 
nitudes of the asymptotic fields are unequal, B^P — 1.5 and 
B% =0.5, where magnetosphere and magnetosheath quan- 
tities are denoted with the superscripts "sp" and "sh", re- 
spectively. The guide field B z is zero. Pressure balance then 
fixes the density profile as a gradient from 1 to 2/3 across 
the current layer with a relative maxima where B x = 0. 

The general features of Figure 1 are consistent with ear- 
lier calculations. One consequence of the asymmetry in field 
strengths across the magnetopause is the evident tendency 
for the magnetic island to grow preferentially in the direction 
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Figure 1. Reconnection in a magnetopause-like configura- 
tion with no guide field. In each panel t = 13.50.^. (a) 
The out-of-plane current density J z , (b) the in-plane field 
lines, and (c) the ion velocity vectors, averaged over a scale 
of ~ di. Note the flows associated with the X-line. The 
longest vector corresponds to a velocity of ~ 0.15tM- 



Figure 2. Cuts at x — 20d; of the outflow region of Figure 
1. Each plot has been averaged over 0.25di in the horizontal 
direction to reduce the noise, (a)-(c) The reconnecting field, 
the out-of-plane field, and the density, respectively. 
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of the magnetosheath. Since equal amounts of flux must re- 



connect from either side in a given time, B a 



B s 



^> v^. The larger inflow on the 



and so vf% = (B sp /B 
magnetosheath side of the reversal surface is evident in the 
ion velocity vectors shown in Figure lc. With the stronger 
inflow, the (frozen-in) magnetosheath field is more strongly 
advected toward the X-line than its magnetospheric coun- 
terpart, leading to the wider opening angle (i.e., the bulge) 
on the magnetosheath side of the current layer. In earlier 
MHD simulations it was shown that the asymmetry of the 
island did not appear when the amplitudes of the reversed 
field on either side of the magnetopause are equal [Scholer, 
1989]. 

Because of the simulation's relatively small size clearly 
separated MHD and kinetic discontinuities do not exist, 
making comparisons with previous work difficult. Neverthe- 
less, in cuts of the reconnection outflow (see Figure 2), the 
self- generated, out-of-plane magnetic field is evident and the 
positive (negative) value of B z on the magnetosheath (mag- 
netosphere) side of the current layer is consistent with earlier 
predictions. The perturbation, however, is much stronger 
on the magnetosheath side; the anti-symmetry of B z seen 
in systems with a symmetric pressure does not extend to 
reconnection at the magnetopause. 

The sharp density drop at the magnetosphere edge of the 
current layer (y w 3.4di) is a consequence of the different 
inflow speeds on the two sides of the layer and the plasma 
mixing across the island. Since the magnetospheric field is 
much stronger, and therefore the inflow velocity from this 
side of the current layer is weak, only a small amount of 
magnetospheric plasma crosses the separatrices. The result 
is an outflow region almost completely comprising magne- 
tosheath plasma. This effect was evident in earlier hybrid 
simulations [Kmuss-Varban et al. 1999; Nakamura and Sc- 
holer, 2000] and has also been seen in satellite crossings of 
the magnetopause [Eastman et al, 1996]. 

3.2. Component Reconnection and Diamagnetic 
Propagation 

In the presence of a magnetic field any nonparallel pres- 
sure gradient produces a diamagnetic drift, 



qjnB 2 ' 



(1) 



where pj = nTj is the thermal pressure and qj is the charge 
of species j. Due to the charge dependence, ions and elec- 
trons drift in opposite directions. To have diamagnetic flows 
parallel to the current layer in our geometry the guide field 
B z and the gradient of the pressure must be nonzero where 
B x =0. 

Even though diamagnetic drifts do not correspond to 
actual particle motions they nevertheless advect the mag- 
netic field [Coppi, 1965; Scott and Hassam, 1987]. In 
two dimensions the magnetic field can be written as B = 
z X Vip(x,y) + B z (x, y)z where ip is the magnetic flux func- 
tion. Taking the cross product of Faraday's law with z yields 
dt'Vip ~ cVS z = 0, or E z = dtip/c. Next, dotting the elec- 
tron fluid equation with z gives 



E z = z • 

c 



X B) - 



m e dv ez 
e dt 



(2) 



The last term represents the inertial effects of electrons and 
breaks the frozen-in condition. Since it is usually small we 
ignore it and substitute for B to get E z — — v e • V^/c and 



a convection equation for the flux [Coppi, 1965], 

d t 1p + V e • W = 0. (3) 

Hence, the electron fluid velocity, which includes a diamag- 
netic component given by (1), advects magnetic structures. 

Strictly speaking, equations (1) and (3) are fluid results 
and need not describe the dynamics at the X-line of our 
simulations where small-scale structures leave the fluid ap- 
proximation on unsure ground. However, as can be seen 
in Figure 3, they are still good descriptions of the system. 
The out-of-plane current density essentially maps the mag- 
netic field lines (see Figure lb) so the island (centered at 
x/di ~ 5) grows robustly during the time shown. Simulta- 
neously, the X-line propagates to the left, the direction of 
the electron diamagnetic drift. The locations of the X-line 
in the three panels of Figure 3 are x/di — 19.2, 14.8 and 
10, respectively, corresponding to a drift speed of O.GIva 
and consistent with the value calculated from the initial pa- 
rameters of v* e — 0.56. As flux reconnects the magnetic 
island widens, decreasing the pressure gradient across its 
center, slowing its drift, and allowing the X-line to overtake 
it. This stabilizes the island's growth because the plasma 
outflow from the X-line towards the near end of the island 
must slow due to the increase in magnetic pressure. The 
actual collision of the X-line with the island, which occurs 
later in time, may be unphysical since in any real system the 
initial separation of the X- and O-lines will be much greater 
than the 12.8di of the present model. Thus, the effects of 
the collision will not be discussed further. 

A surprise is that the diamagnetic propagation evident in 
Figure 3 has not been seen in previous hybrid simulations 
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Figure 3. The out-of-plane current density at three times 
for a magnetopause geometry similar to Figure 1, but in- 
cluding an out-of-plane guide field and therefore an in-plane 
diamagnetic drift. The initial conditions are those described 
in section 2.2 with an asymptotic guide field of 3.0 in the 
magnetosheath and 3.8 in the magnetosphere. The density 
drops from 1 to 0.1 across the current layer. 
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of the magnetopause [Krauss-Varban et ai, 1999; Nakamura 
and Scholer, 2000] even though these models should include 
the relevant electron diamagnetic drifts. In some cases the 
simulation parameters may be such as to minimize diamag- 
netic effects {e.g. a lower magnetospheric temperature). 
Diamagnetic drifts may also be absent if the simulations 
were performed with an initially uniform out-of-plane field. 
According to the force balance condition, the local pressure 
gradient (and hence the drift) must then be zero at the re- 
versal surface, although whether it would remain zero as 
the system evolves is unclear. In section 5 we discuss fur- 
ther why a uniform guide field may not be representative 
of reconnection at the magnetopause. A final possibility 
occurs for simulations using a fixed, spatially- localized re- 
sistivity to break the frozen-in condition at the X-line. How 
such a resistivity model would affect an X-line that is trying 
to propagate is unclear, but it is nevertheless evident from 
our results that such models may be inappropriate for sim- 
ulating guide field reconnection at the magnetopause when 
diamagnetic drifts are large. 

Figure 4a shows the current density for the simulation of 
Figure 3. Unlike the example shown in Figure 1, the mag- 
netic island does not bulge into the magnetosphere side of 
the current layer. This is because the reconnecting magnetic 
field B x is anti-symmetric across the current layer and thus, 
as discussed earlier, both sides contribute equal amounts of 
magnetic flux to the reconnection. Clearly evident in Figure 
4a, as well as the last two panels of Figure 3, is the left-right 
asymmetry of the opening angle of the magnetic field lines 
near the X-line. This is caused by the interplay of the X- 
line's diamagnetic drift, controlled by just the electrons, and 
the reconnection outflow, which must also involve the ions. 

Consider the relative motion of the X-line and the ions. 
Downstream from the reconnection site the J X B force ac- 
celerates frozen-in ions up to the Alfven speed, iiVAx- To 
this must be added the relative velocity between the diamag- 
netic drift of the ions and the X-line, u* = \v*i\ + \v* e \, with 



^ 6 



4 



10 



6 



4 





1 1 1 1 


1 1 i 1 


1 1 1 


i i i i 


1 1 u 

lal 







5 


10 


15 

x/d, 


20 


25 












1) 













































5 


10 


15 


20 


25 



x/dj 



Figure 4. Structure of (a) the out-of-plane current density 
J z and (b) ion flow vectors at t = 16.0f2~ for the simulation 
shown in Figure 3. The large flows at the current layer are 
the ion diamagnetic drifts. The longest vector corresponds 
to a velocity of ~ 1.2va 



the result being a right-left asymmetry in the outflow veloc- 
ities, ivAx + v*. Since, from continuity, the opening angle 
of the magnetic field downstream from the X-line is roughly 
Vin/v ou t — and Vi„ is roughly constant due to the symmet- 
ric reconnecting field — we expect the rightward opening 
angle to be much smaller than its counterpart, consistent 
with Figure 4b. 

By a similar argument, normal ion outflow from an X- 
line only develops when vax > «*■ For the simulations pre- 
sented in Figures 3 and 4 the electron and ion diamagnetic 
drifts at the reversal surface are —0.56 and 1.12, respec- 
tively, compared with vax = 1-35 (based on a density of 
0.55, corresponding to a mixture of the magnetopause and 
magnetosheath plasma). Thus, for the parameters of this 
simulation, the reconnection should be marginally unable to 
reverse the ion flow to the left of the X-line. The ion flow 
vectors are shown in Figure 4b at the same time as the cur- 
rent density shown in Figure 4a. The strong ion flow on 
the right side of the X-line results from the combination of 
the ion diamagnetic velocity with the direct acceleration by 
the reconnected magnetic field. The reconnection generated 
forces almost reverse the diamagnetic flow just to the left of 
the X-line. 

We again emphasize that it is the relative drift of the 
ions and electrons, i.e. v, = \v*i\ + \v* e \, that must be 
smaller than the in-plane Alfven speed for normal reconnec- 
tion to proceed. Reconnection was suppressed in a nearly 
indistinguishable fashion in a series of simulations with the 
same total temperature T e + Ti (and therefore identical v*) 
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Figure 5. Cuts at x = 13.8di of the outflow region of Figure 
4. Each plot has been averaged over O.ldi in the horizontal 
direction to reduce the noise, (a)-(c) The reconnecting field, 
out-of-plane field and density, respectively. 
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but differing values of T e and T t (Ti = 1/2, T e = 5/2 and 
Ti = 5/2, T e = 1/2). Since T » T e at the magnetopause, 
stabilization can still occur even if the drift of the X-line, 
which is due only to v te , is negligible. 

The cuts in Figure 5 show that for this simulation the re- 
connection field B x is symmetric across the layer while the 
pressure drop across the magnetopause is supported by the 
jump in the guide field B z . In contrast to the density pro- 
file in the case with no guide field, the density drop across 
the layer is nearly centered on the magnetic island. Large 
drops occur on either edge of the current layer with a weak 
plateau in the middle of the island. The inflow velocity into 
the island from either side of the current layer is now roughly 
equal because of the near anti-symmetry of B x . The plasma 
in the core of the island is therefore a nearly equal mix of 
magnetosphere and magnetosheath plasma and is not dom- 
inated by the magnetosheath, in contrast to the case shown 
in section 3.1. 

4. Suppression of Component Reconnection 

We have shown for a single case that, consistent with 
equation (3), X-lines convect with the diamagnetic velocity 



of the electrons in the reversal region and suggested a cri- 
terion for the suppression of reconnection from diamagnetic 
drifts, v, > vax- We now more fully explore this stabiliza- 
tion criterion by varying the strength of the out-of-plane 
field to vary the ratio of w» to vax- For our magnetopause 
model the diamagnetic velocity at the surface where B x — 
can be written as 



cT, 



qjLpjBz 



(4) 



where the pressure scale length is given by L~ 1 = \d y p\/p. 
Varying the amplitude of the out-of-plane field B z while 
keeping all the other parameter fixed alters v*j but keeps 
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Figure 6. The top panel shows the position of the X-line 
versus time for several runs differing only in the strength of 
the guide field B z . The solid line is for the run of section 
3.1 and experiences no diamagnetic motion. The dotted, 
dashed, and dash-dotted lines correspond to B z q = 3, 1.5, 
and 0.5, respectively. The bottom panel is a comparison of 
the measured X-line velocities and the electron diamagnetic 
velocity at the reversal surface. The straight line has a slope 
of one and is plotted to guide the eye. In both panels the 
X-point location corresponds to a saddle point in ip. 
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Figure 7. Reconnected flux versus time for four simula- 
tions. The reconnection rate is the slope of each curve. The 
solid line is a reference run with no diamagnetic drift. The 
dotted, dashed, and dash-dotted lines are from simulations 
with differing values of B z o (1-5, 1.0 and 0, respectively) and 
therefore v*. 



4 



10 



>^ 6 
4 



- 1 1 












1 1 u 











5 


10 


x/d, 


15 


20 


25 














1) 













































































5 


10 


x/d, 


15 


20 


25 



Figure 8. (a) The out-of-plane current density J z and (b) 
ion flow vectors at t = 13.0f2~ from a simulation with an 
asymptotic B z = 0.5 and a strong diamagnetic drift. The 
large flows to the right at the current layer are the result of 
the ion diamagnetic drift. The longest vector corresponds 
to a velocity of ~ 0.95ua 
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vax constant. As can be seen in Figure 6a, a large guide 
field implies a small v* . Despite the occasional glitches (dis- 
continuities in the motion arise when a new X-line triggered 
by random perturbations begins to dominate), the drift ve- 
locity is remarkably constant. 

The bottom panel of Figure 6 shows the initial elec- 
tron diamagnetic velocity at the reversal surface versus the 
drift speed of the X-line measured in the simulation. Al- 
though the agreement for sub-Alfvenic velocities is quite 
good, as v te approaches vax another effect becomes impor- 
tant. Rather than move the ions at super- Alfvenic velocities 
the system develops an electric field E y transverse to the cur- 
rent layer that, combined with the guide field B z , generates 
an E X B drift opposing the ion motion and adding to the 
electron velocity. There is no change in the net current but 
the electrons, and hence the X-line, move faster than just 
the diamagnetic velocity would suggest. 

Since diamagnetic flows of sufficient strength disrupt the 
large-scale flows characteristic of reconnection, increases in 
v, should correlate with decreases in the amount of recon- 
nected flux. Figure 7 depicts the extent to which a diamag- 
netic drift can alter the reconnection rate of a system. The 
reconnected flux is plotted versus time for four different sim- 
ulations: one with no diamagnetic drift, as in Figure 1, and 
three with varying asymptotic guide fields, B z o = 1.5, 1.0, 
and 0. For the latter cases the associated drift speeds (dia- 
magnetic plus E X B) at the reversal surface are v, = 0.8, 1.1 
and 1.5, respectively. Reconnection is almost completely 
suppressed for the largest drifts. The plateau in reconnected 
flux seen in each of the bottom three curves is an artifact 
(discussed earlier) of our system's periodic boundary condi- 
tions that allow the X-line to impinge on the more slowly 
propagating magnetic islands. 

For large guide fields the diamagnetic drift speed ap- 
proaches zero and the stabilization is minimized. The con- 
verse limit is more subtle since pressure balance imposes 




Figure 9. Schematic of the asymptotic magnetic fields 
in the magnetosphere, B sp , and magnetosheath, B sh . The 
reconnecting and out-of-plane components of the magnetic 
field depend on the orientation of the reconnection plane: 
x — z versus x' — z'. 



restrictions on the system's configuration. If the profiles 
of the density and the reconnecting field are specified then 
there is a minimal possible value of the guide field in the 
current layer; B z cannot be taken to be zero and the drift 
cannot grow without bound. Hence although (for example) 
the GEM challenge simulations had zero guide field they 
did not have a large diamagnetic drift. Rather, B z — > and 
L p — > oo, keeping v* = 0. 

To confirm the hypothesis that the reconnection process 
cannot reverse the ion diamagnetic drift and drive outward 
flows from the X-line for large values of the diamagnetic ve- 
locity, we show J z and the ion velocities from a simulation 
with an asymptotic out-of-plane field B z — 0.5 in Figure 
8. In this case the island of reconnected flux remains small 
even at late time because of the weak growth shown in Fig- 
ure 7. The flows are actually smaller than the peak of the 
expected diamagnetic velocity because of the development 
of the electric field E y across the layer. A slight reduction 
of the ion flow speed to the left of the X-line indicates that 
the reconnected field lines exert a small leftward directed 
force on the ions in this region, but the result is only a small 
perturbation of the ion flow pattern. Thus, the X-line is 
unable to couple effectively to the ions, strong reconnection 
does not develop, and the island amplitude saturates. 

5. Dependence on the Reconnection Plane 

Although two-dimensional simulations are computation- 
ally cheaper than their three-dimensional counterparts, this 
simplicity comes with drawbacks. Perhaps the most impor- 
tant when considering diamagnetic drifts is that the orien- 
tation of the X-line with respect to the magnetic field con- 
figuration is externally imposed, rather than being free to 
develop. Different simulation planes have different recon- 
necting and out-of-plane fields and hence varying amounts 
of diamagnetic stabilization. In this section we explore this 
effect, albeit only by considering rotations around the y-axis 
to avoid introducing stabilizing B y components into the sys- 
tem. We will assume that the favored plane (i.e., the one 
the system would choose if allowed to evolve in three dimen- 
sions) is the one where the reconnection rate is maximal, but 
note that the results are not a replacement for more ambi- 
tious 3-D simulations. We further assume that the diamag- 
netic effects are playing the dominant role in determining 
the plane where reconnection is most robust. Thus, we are 
considering systems in which v* ~ vax- 

Three-dimensional particle simulations have already ad- 
dressed this issue for the basic case of a system with no 
guide field and no pressure drop across the current layer 
[Hesse et ai, 2001; Pritchett, 2001]. Random perturbations 
were used as seeds to avoid imposing a plane a priori. Al- 
though X-lines were free to develop with any orientation, 
they preferentially formed parallel to the current layer, con- 
sistent with our assumption of treating only rotations about 
y. Moreover they remained quasi- two-dimensional even at 
late times. However, the simplicity of B makes this case de- 
generate since the reversal surface (the place where B x — 0) 
must be at the center of the current layer for any rotation 
about y. This restriction does not exist in more complicated 
configurations. 

The basic geometry of magnetic fields at the magne- 
topause is shown in Figure 9. In order to agree with our sim- 
ulation coordinates, the magnetospheric field B sp is taken to 



8 



SWISDAK ET AL.: DIAMAGNETIC SUPPRESSION OF RECONNECTION 



be in the positive x direction while the magnetosheath field 
B is inclined at an angle 8 with respect to the negative 
x axis. Thus, the rotation angle across the magnetopause 
is 7r — 6. The y direction is perpendicular to the surface of 
the magnetopause. For a two-dimensional simulation in the 
x — y plane the reconnecting components of the magnetic 
field will be parallel to x, 



g s v _ B "p 

B sh = _ B sh CQse 



and the guide field to z, 



B sh sml 



(5) 
(6) 



(7) 
(8) 



The primed axes in Figure 9 are rotated clockwise by an 
angle a with respect to the unprimed system and represent 
another possible reconnection plane. In a simulation carried 
out in the x' — y' plane the reconnecting components of the 
magnetic field will be 



b sp = B °P cosa 
B sh = _ B °h cos (0 + a ^ 

while the guide field components are 
B sp = _ B s P ginQ 

B$ = B sh sm{6 + a). 



(9) 
(10) 



(11) 
(12) 



As shown in section 4, stabilization via a diamagnetic 
drift is inversely proportional to the value of the guide field 
B z in the center of the current layer. Thus, the expectation 
is that the plane with the maximum value of the guide field 
at the surface where reconnection takes place will dominate. 
Unfortunately, the asymptotic fields determine neither the 
location of the reversal surface nor the value of the guide 
field there. Furthermore, the actual guide field at the rever- 
sal surface may vary as reconnection proceeds. To simplify 
matters, we adopt the crude hypothesis that the guide field 
at the reversal surface B g is the mean of the asymptotic 
guide fields, which is [B sh sin(# + a) — B sp sina]/2 in the 
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Figure 10. Reconnected flux versus time for three simula- 
tions. The solid line is a simulation with an asymptotic guide 
field of 1.5 (the dotted line in Figure 7). The dashed and 
dotted lines are simulations in planes rotated ±15° about y. 



rotated system. With this ansatz the optimal plane for re- 
connection is that where B g is maximized with respect to 
the angle a, from which we derive the condition 

B sp cosa = B 3h cos(6 + a). (13) 

Interestingly, this corresponds to the plane where the recon- 
necting field components are equal in magnitude but oppo- 
sitely directed (see equations (9)-(10)). Figure 10 shows the 
effect of rotating our simulation plane on the reconnection 
rate and suggests that when diamagnetic drifts are strong, 
reconnection in the anti-symmetric plane is favored. 

6. Summary and Discussion 

Diamagnetic drifts can be an important part of mag- 
netopause reconnection in certain parameter regimes. In 
our simulations the plasma (5 (approximately 1) and the 
magnetosheath-magnetosphere density contrast (a factor of 
10) are roughly consistent with observational constraints. 
However our initial conditions also neglect both tempera- 
ture gradients and scale lengths larger than an ion inertial 
length. For other magnetopause models that include these 
features the diamagnetic drifts and associated stabilization 
may be smaller. 

When diamagnetic drifts are significant an X-line will ad- 
vect with the electron fluid velocity, causing the separatri- 
ces to have an asymmetry in their opening angles with the 
wider facing the direction of motion. For large enough flows, 
v* — \v* e \ + \v*i\ ~ va, reconnection can be completely sup- 
pressed since the large-scale outflows from the X-line needed 
for fast reconnection cannot develop. The reconnected field 
lines are not strong enough to reverse the ion diamagnetic 
flows toward the X-line. 

The diamagnetic effects presented here introduce a (3 de- 
pendence that was absent in the symmetric systems studied 
earlier. The diamagnetic stabilization condition v* > va 
can be rewritten as a condition on /3: 



> 



B x c/uj p 



(14) 



where B x is the amplitude of the reconnecting field, fi x — 
8wnT/ B x , and B z and L p are evaluated at the current layer. 
For B z ~ B x and for pressure scale lengths of the order of 
the ion inertial length c/to p i, the threshold for the suppres- 
sion of reconnection is (3 ~ 1. 

A statistical study of accelerated flow events was done by 
Scurry et al. [1994] to investigate the dependence of day- 
side reconnection on magnetosheath j3 and the shear angle 
of the magnetic field across the magnetopause. They found 
that for low values of (3 the accelerated flow events spread 
over a broad range of clock angles while at high values of f3 
they were strongly correlated with large clock angles. The 
implication is that magnetic reconnection in the presence of 
a significant guide field is suppressed at high /3. Thus, these 
observations are consistent with the criterion given by equa- 
tion (14) for the diamagnetic suppression of reconnection. 

When B z 3> B x equation (14) indicates that it is difficult 
to stabilize reconnection with diamagnetic effects, just as is 
demonstrated in section 4. We do not address what happens 
when the guide field becomes very large, but Rogers et al. 
[2001] examined the effect of /3 on reconnection in a sym- 
metric system (no pressure drop across the current layer) us- 
ing theoretical arguments as well as numerical simulations. 
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They argued that fast reconnection as described in section 
3 can take place whenever dispersive waves (e.g., whistlers 
or kinetic Alfven waves) exist in a system at small spatial 
scales. Only under extreme conditions (B z ~ B x m; /m e , 
(3 X = 8imT/B% <JC 1) is reconnection inhibited by a guide 
field B z . 

A final curious feature of equation (14) occurs in the limit 
where B z is small, as the expression would seem to imply 
that reconnection is inhibited even for low values of j3 x . This 
is not the case because there is a linkage between the pres- 
sure scale length L p and the value of B z at the reversal 
surface. Since the reconnection field B x = at this lo- 
cation, pressure balance requires that p/L p = B Z B' Z with 
B' z — \d y B z \. This condition implies that B Z L P = p/B' z , 
and so that if the right-hand side is to remain finite then 
L p — » oo as B z — ► 0. In the absence of a guide field the 
pressure gradient must be zero at the reversal surface. For 
nonconstant B z , the minimal initial value of the reversal 
guide field is determined by the details of the initial condi- 
tions. 

Unlike systems with simple reversed fields, reconnection 
at the magnetopause is not limited to a single plane, a fact 
that has led to the prediction that the magnetopause field 
might become stochastic [Galeev et at, 1986; Lee et al, 
1993]. Diamagnetic drifts offer a way to trigger reconnec- 
tion at multiple surfaces since they can directly affect which 
planes dominate reconnection. A rough argument suggests 
that the maximal guide field (and hence the minimal drift 
and minimal stabilization) occurs at the surface where the 
components of the field that reconnect are equal and oppo- 
site. Hence we expect this to be the dominant location of 
"component reconnection" at the magnetopause. 

The direction of the X-line advection can in principle be 
predicted from the local pressure gradients and magnetic 
fields. Under appropriate conditions this effect should be de- 
tectable by spacecraft at the magnetopause, although defini- 
tive results are certainly more likely with multiple spacecraft 
missions like Cluster or the future Magnetospheric Multi- 
scale Mission. 
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